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During a component compatibility test of the NASA HiVHAc Hall thruster, a number of 
plasma diagnostics were implemented to study the effect of varying facility background 
pressure on thruster operation. These diagnostics characterized the thruster performance, 
the plume, and the plasma oscillations in the thruster. Thruster performance and plume 
characteristics as functions of background pressure were previously published. This paper 
focuses on changes in the plasma oscillation characteristics with changing background 
pressure. The diagnostics used to study plasma oscillations include a high-speed camera and 
a set of high-speed Langmuir probes. The results show a rise in the oscillation frequency of 
the “breathing” mode with rising background pressure, which is hypothesized to be due to a 
shortening acceleration/ionization zone. An attempt is made to apply a simplified ingestion 
model to the data. The combined results are used to estimate the maximum acceptable 
background pressure for performance and wear testing. 
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Glenn Research Center 
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= Mean Channel Diameter 
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Science Mission Directorate 
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= Charge Exchange 
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In-Space Propulsion Technology 

SEE 

= Secondary Electron Emission 

HiVHAc = 
EDU 

High Voltage Hall Accelerator 
Engineering Development Unit 

Cl 

= Confidence Interval 


I. Introduction 

T HIS paper is a continuation of a previous work on the effect of background pressure on the performance and 
plume of the NASA High-Voltage Hall Accelerator (HiVHAc) thruster. 1 While the previous paper focused on 
the breakdown of performance into various efficiency factors and far-field time-averaged diagnostics, this paper will 
focus on plasma oscillation characteristics and time -resolved diagnostics. 

The HiVHAc project is funded by the NASA Science Mission Directorate’s (SMD’s) In-Space Propulsion 
Technology (ISPT) program with the goal of delivering a Hall thruster system beyond the current state-of-the-art 
(SotA). The project objective is to enable various Discovery class missions while keeping the cost of the propulsion 
system low relative to other electric propulsion systems. Mission studies performed in 2009 showed that a 3.5-kW 
HiVHAc system is able to outperform an SotA 4.5-kW flight Hall thruster by delivering 6-12% more payload mass 
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on four missions of interest. 2 These missions include Vesta-Ceres rendezvous (the Dawn mission), Koppf comet 
rendezvous, Nereus (a near-Earth asteroid) sample return, and NEARER, which involves two near-Earth asteroid 
returns. The thrust profiles of the missions range from ones that favor high specific impulse to ones that favor high 
thrust-to-power demonstrating the wide throttle -ability of the 3.5-kW HiVHAc. 

The iteration of the HiVHAc Hall thruster tested in this work is the Engineering Development Unit 2 (EDU2), 
which is designed to operate at up to 3.9 kW discharge power. This version of the HiVHAc thruster features a 
discharge channel replacement mechanism designed to achieve long operational life time for high specific impulse 
(high discharge voltage) operation. Compared to the previous iteration, which is the EDU1, the EDU2 has better 
thermal management, a simpler design for the channel replacement mechanism, and improved voltage isolation. The 
EDU2 has undergone the performance acceptance test (PAT), 3 a vibration test, 3 a plume divergence characterization 
test, 4 in-situ mechanism actuation tests (unpublished), and a component compatibility test (CCT). 5 The background 
pressure study described in this paper was performed during the CCT in April and May of 20 13. 5 This work 
produced data that can guide future HiVHAc testing, including the wear test. 

To operate in a space-representative vacuum environment, the EDU2 was tested in NASA Glenn Research 
Center’s (GRC’s) Vacuum Facility 5 (VF5). Many diagnostics were deployed as a part of the testing, including a 
thrust stand, a retarding potential analyzer, a Wien filter, also known as an ExB probe, a Faraday probe, an infrared 
camera, a high-speed camera, and a set of high-speed Langmuir probes. Additional details regarding the CCT setup 
can be found in a publication by Kamhawi. 5 The present paper will focus on setup and analysis of the high-speed 
camera and the high-speed Langmuir probes, as well as the implications of the background pressure study for future 
testing. 


II. Past Work 

Many kinds of plasma oscillations have been known to exist in Hall thrusters. 6 Since the 1960’s, many studies 
have been performed to characterize these oscillations. 7 " 20 Yet the roles of these oscillations in relation to thruster 
performance and life-time are still largely unknown. Currently, there are two oscillation modes that are believed to 
be important to Hall thruster operations, the “breathing” mode and the “spokes” mode. Both modes are believed to 
be tied to ionization processes in the Hall thruster because their frequencies are in the kilohertz to tens of kilohertz 
range, which is most closely related to the characteristic ion and neutral transit time across the discharge axial length 
scale. 6, 12 Recently, there has been increased interest in these plasma oscillations because they have been observed 
to vary greatly with facility conditions 18 and may be responsible for differences in thruster behavior between ground 
tests and on-orbit. As more Hall thrusters are being considered for flight applications, including the HiVHAc series, 
it has become critical for the Hall thruster community to better understand the behavior of these plasma oscillations 
so that on-orbit performance can be predicted using ground test data. For the present study, the high-speed camera 
and Langmuir probe were chosen to study these plasma oscillations. 

The high-speed camera has only recently been used as a quantitative tool for studying Hall thruster plasma 
oscillations. High-speed optical imaging on a Hall thruster was first performed by Darnon, et. al., though the view 
was not perpendicular to the annulus plane and did not capture azimuthal plasma movement. The instrument 
Darnon, et. al., used was a charge -coupled device with an image intensifier. 1 1 Subsequently, Lobbia, Liu, and 
Gallimore studied a cluster of annular Hall thrusters using a commercial high-speed camera 21 while Parker, Raitses, 
and Fisch studied a cylindrical Hall thruster. 22 Shortly after, McDonald developed a Fourier-transform technique to 
study both breathing and spokes mode in annular Hall thrusters. 19 Sekerak then combined the Fourier-transform 
technique with results from high-speed Langmuir probes to study oscillation mode transition in Hall thrusters as a 
function of magnetic field strength. 20 The analysis technique used in the present study is based on the 
aforementioned Fourier-transform technique. 

The high-speed Langmuir probe was first used on a Hall thruster by Lobbia and Gallimore, who first identified 
all of the distortion effects associated with high-speed sweeping of Langmuir probes up to the megahertz range. 23 
Lobbia then combined a high-speed dual Langmuir probe (HDLP) with a Fourier-transform correlation technique to 
temporally and spatially resolve the breathing mode in a Hall thruster. 17 The high-speed Langmuir probe technique 
used by Sekerak is a refined version of Lobbia’ s work. 20 

The breathing mode oscillation is believed to be an axial ionization wave. The name comes from the 
resemblance of the exhausted plasma packets to that of breaths of gas during breathing mode oscillation. The 
simplest quantitative description of the breathing mode oscillation is based on modeling work done by Fife. Fife 
started with the assumption that the breathing mode is an axial ionization-driven oscillation and derived the scaling 
relationship shown in Eq. (1) below. 12 
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In Eq. (1), fi is the frequency of the breathing mode, u* is the speed of the ions leaving the thruster, u n is the speed of 
the neutrals entering the ionization zone, and Li is the characteristic length over which the oscillation action is taking 
place. Typically, Li is taken to be the overall length of the discharge zone, where both ionization and acceleration are 
taking place. This scaling relationship has been shown to predict general trends in the breathing mode oscillation 
frequency with discharge voltage though the exact values do not always agree. 12 The role of the breathing mode in 
electron transport is currently unclear, though Lobbia has argued that the breathing mode is an intrinsic feature in the 
Hall thruster and is necessary for the thruster to achieve enough electron transport to maintain the discharge. 17 The 
scaling relation was not created with the ingestion of background neutrals in mind. 

The spokes mode oscillation is believed to be an azimuthal ionization wave. The name comes from the 
appearance of rotating plasma packets, which resemble spokes on a wheel, in the channel of a Hall thruster when 
viewed from downstream of the exit plane. Works performed by Janes and Lowder 7 and McDonald. 19 have provided 
some basic understanding of the spokes mode oscillation. The spokes mode is presented as an alternate means by 
which electron current can be transported across magnetic field lines. Subsequent work by Sekerak showed that a 
Hall thruster can operate with almost no breathing mode and strong spokes mode oscillations present. Sekerak 
observed that the presence of spokes mode without breathing mode is correlated with high thrust efficiency for a 
medium-power Hall thruster. 20 Spokes have been observed to travel in the same direction as the Hall current and 
anywhere from two to six spokes can be observed simultaneously. 19, 20 There are currently no existing scaling 
relationships for the number and speed of the spokes present. 

III. Experimental Setup 

The experimental setup was previously reported 1 and is repeated here with new details about the high-speed 
diagnostics for convenience. To simplify plot labeling, operating conditions will be labeled as vvv-k.k-Pnx, where 
vvv is the discharge voltage in volts, k.k is the discharge power in kilowatts, and n is the normalized background 
pressure, defined as the number of multiples of the lowest achievable operating background pressure for the throttle 
point vvv-k.k. This pressure is measured by an ion gauge located directly below the thrust stand. Lor example Plx 
indicates the operating condition corresponding to the minimum achieved background pressure, while PI Ox 
indicates the background pressure is ten times that of Plx. Unless otherwise noted, all spatial positions presented in 
this paper have been normalized by the mean channel diameter (MCD) of the thruster. MCD is defined as the 
average of the inner and outer discharge -channel wall diameters. 

A. Test Article and Test Matrix 

The NASA HiVHAc EDU2 is a 3.9-kW xenon Hall 
thruster. The thruster is highly throttle -able. The discharge 
voltage can vary from 200 to 650 V, corresponding to specific 
impulse values of 1200 to 2700 s. The discharge power ranges 
from 0.3 to 3.9 kW. During testing, the magnetic field settings 
are selected to create a magnetic lens symmetric about the 
channel centerline, while maximizing total efficiency. This 
magnetic field setting optimization is performed at the lowest 
achievable background pressure for each throttle point, figure 1 
shows a photograph of the test equipment around the NASA 
HiVHAc EDU2 thruster. 

Lor the purpose of this test, seven throttling points were 
selected. These seven points are 300-1.5, 300-3.0, 400-2.0, 400- 
3.2, 500-2.0, 500-3.0, and 500-3.9. At each throttling point, the 
thruster settings (including magnet currents) are optimized at 
the lowest achievable pressure as measured by the ion gauges 
next to the thrust stand. Then, the background pressure was 
raised by injecting xenon via an auxiliary flow line that exits at 
>4 meters downstream of the thruster pointed away from the 
thruster. 



Figure 1. A photograph of the test equipment 
around the NASA HiVHAc EDU2 thruster. 


3 

American Institute of Aeronautics and Astronautics 




B. Test Facility 

Testing was conducted in NASA GRC’s VF5. This cylindrical facility is 4.6 m in diameter, 18.3 m long, and is 
pumped with cryo-panels and 20 oil diffusion pumps. The thruster was mounted on a thrust stand located close to 
the cryo-panels, with the thruster firing away from the panels. Facility pressures were monitored with four ion 
gauges, three of which were mounted next to the thrust stand and the fourth mounted on the facility wall at the same 
axial location as the thruster. Manufacturer specifications state that the ion gauges are accurate to ±6% of reading. 
The positions of the gauges are shown in Fig. 1. Ion gauges 1 and 2 are both facing downstream while ion gauge 3 is 
facing upstream. Ion gauge 1 and 2 agree to within 10% of each other. Ion gauge 3 reports 0.63 to 0.72 times the 
reading of ion gauge 2. Ion gauge 2 readings were used to determine the number of multiples of the lowest 
achievable background pressure that the thruster was experiencing. The lowest pressure achieved for each tested 
condition, corrected for xenon, varied from 1.4xl0" 6 to 2.2xl0" 6 Torr. 

Thruster propellant was supplied via commercially available mass flow controllers. These mass flow controllers 
were calibrated using xenon prior to testing. Typical uncertainty of measurement was ±1.0% of reading. Research- 
grade xenon was used. 

Thruster electrical power was supplied via the Colorado Power Electronics (CPE) brassboard #2 (BB#2) power 
processing unit (PPU). 24 This PPU is capable of supplying 200 to 700 V at up to 4 kW over the main discharge line 
as well as supply power to the cathode heater, keeper, and magnet electrical circuits. The PPU was located outside 
of the vacuum facility during testing. 


C. High-Speed Diagnostics 

The high-speed camera was a commercial device on loan from the Aerospace Corporation. To obtain a non- 
oblique view of the thruster channel without putting the high-speed camera into the plume a large plane mirror was 
deployed. The camera was positioned outside of the vacuum chamber and viewed the mirror through a window port. 
The camera was set to capture at a frame rate of roughly 80,000 frames per second at 256 by 256 pixels. In order to 
capture the cathode as well as compensate for possible misalignment during pump down the view area was larger 
than the discharge channel area, which took up a 100-by- 100-pixels area. Figure 2 shows an overhead diagram of 
the positions of the mirror and camera. The view direction of the camera was actually 4° off-axis rather than along 
firing axis due to the need to accommodate an IR camera. This small offset was considered acceptable as it produces 
negligible impact on the data analysis. High-speed video at varying background pressure was only collected for the 
throttle points 500-2.0, 500-3.0, and 500-3.9. High-speed video was obtained at only the lowest achievable pressure 
for the throttle points 300-1.5, 400-2.0, and 400-3.2. 

The high-speed Langmuir probe array consisted of four pairs of Langmuir probes. Each pair consisted of one 
active probe and one null probe. Data from the null probe was used to correct for capacitance effects, which allows 
for accurate high-speed Langmuir probe sweeps. The entire probe array was mounted on a two-axis motion system 
and was well outside the plume when not in use. The four pairs of probes were divided into two groups. The 
upstream group consisted of two pairs set 2 cm apart in the azimuthal direction. The downstream group consisted of 
the other two pairs set 6 cm apart in the azimuthal direction. The axial distance between the upstream and 
downstream groups was 10 cm. There were three measurement positions for the array at each thruster condition. For 
these three positions, the upstream probe was 10, 15, and 20 cm away from the thruster exit plane, respectively, and 
the probe tips were radially aligned with the centerline of the channel. To correlate the Langmuir probe signals 
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Figure 2. A diagram of the diagnostics setup. 


Path of the 
downstream group ~~ ^ 

Two pairs in 
each group 





Path of the _f 


— 30 cm 

— 25 cm 

— 20 cm 

— 1 5 cm 

— 1 0 cm 


upstream group 
Top View 




Figure 3. A diagram of the movement 
path of the high-speed Langmuir 
probe array. 
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across different locations, the discharge electrical line was monitored. Figure 3 shows an overhead diagram of the 
movement path of the probe array and the interrogation 
locations with respect to the thruster. 

A specialized high frequency power supply was used to 
supply a 200-kHz signal through custom circuitry to sweep 
the probes. The custom circuitry was tuned to minimize 
capacitive effects associated with high-speed probe 
sweeping. A custom computer with high-speed data 
acquisition systems, high processing power, and high 
storage capacity is used to record the data. Figure 4 shows a 
photograph of the high-speed Langmuir probe array moving 
into place to record data while the thruster was operating. 

IV. Data Analysis 

A. High-Speed Camera Analysis 

For this study, the data analysis focused on the breathing 
and spokes mode oscillations. The high-speed camera 
videos were processed in the following steps. 

1 . The image was cropped so that only the annulus region was analyzed. 

2. The video was averaged over time to create an averaged image. 

3. A circle fit was performed on the averaged image to find the center of the thruster and the boundaries of the 
discharge. 

4. The data was divided into 120 azimuthal bins and the average intensity in each bin on each frame was calculated. 

5. The binned data was normalized against the time-averaged values of the binned data. 

6. A two-dimensional (2D) Fourier transform was performed to obtain a series of power spectra categorized by the 
mode number m. 

The data resulting from step 5 were the relative fluctuation of each azimuthal bin with time. This step produced 
results that are much less sensitive to variations in pixel sensitivity to light, transmission of the viewport, and 
reflectance of the mirror. Any breathing or spokes mode present showed up more clearly when this step was 
performed. Note that step 5 also removed any real, time-averaged, azimuthal variation in light intensity. For 
instance, the presence of a stationary hot spot would be removed by this step. This shortfall was considered 
acceptable because the present study was focused on oscillatory instead of stationary features. For the power spectra 
obtained from step 6, the m = 0 spectrum showed the frequency of any breathing mode present, the m = 1 spectrum 
showed the frequency of any spokes mode with one spoke present, the m = 2 spectrum showed the frequency of any 
spokes mode with two spokes present, and so on. 

To properly interpret the power spectra from the 2D 
Fourier transform, one must first understand what the 
transform was really doing to the data. The transform broke 
down the data into a series of 2D sinusoids that when taken 
together could re-constitute the original data. Since one of 
the dimensions of the binned data was time and the other 
was azimuthal angle, the transform looked for waves that 
were both repeating in time and in angular space. Figure 5 
shows illustrations of the ideal modal patterns that the 
transform was looking for when processing the data. Parts 
of the annulus with higher than average light intensity are 
shown in white while parts with lower than average light 
intensity are shown in black. Although the true ideal pattern 
have an intensity that varies in a sinusoidal fashion as 
opposed to having discrete black and white shades, this 
simplified illustration is used for better clarity. 

An important feature of the modal patterns in Fig.5 is that the angular length of the light packet is necessarily 
shorter for higher mode number patterns than for lower mode number patterns. For example, the light packets of m 

5 

American Institute of Aeronautics and Astronautics 


" = 1 = Q'Q Q> Q 

o o © 

Figure 5. Illustrations of ideal modal patterns as 
applied in the 2D Fourier transform analysis. 



Figure 4. A photograph of the high-speed 
Langmuir probe array with the thruster 
operating. Only the downstream group of probes 
is visible due to the camera angle. 



= 3 mode are 60 degrees long while the light packet of m = 1 mode is 1 80 degree long. Another important feature is 
that these model patterns have light packets that are equally spaced. However, real Hall thruster plasma waves are 
generally more random, with light packets that can be shorter or longer than ideal and with spacing that are not even. 
Thus, a spokes mode oscillation pattern with two short light packets may be de-composed into a series of modal 
patterns at the m = 4 mode in addition to the m = 2 mode because the length of the light packets are closer to that of 
the ideal m = 4 modal pattern but the spacing of the packets are more similar to the ideal m = 2 modal pattern. 
Another example is when 4 light spokes are present with uneven lengths and spacing, and is decomposed into a 
series of peaks at m = 3, 4, and 5. The real spokes pattern did not fit neatly into the ideal modal patterns but that 
does not necessarily mean there were 3, 4, and 5 light spokes present at the same time; instead, the real pattern is a 
super-positioning of m = 3, 4, and 5 modal patterns on top of each other. To complicate matters further, if the 
thruster changes from an ideal m = 3 modal pattern to an ideal m = 4 modal pattern, and then changes to an ideal m 
= 5 modal pattern during the filming of the high-speed video, the resulting Fourier transform can also produce a 
series of peaks at m = 3, 4, and 5. Thus, close examination of the actual high-speed video is sometimes necessary to 
determine the exact nature of the modal pattern present. 

In contrast, analysis of the breathing mode oscillation can be completed by either averaging the light intensity 
data across all bins and performing a ID Fourier transform or taking the m = 0 result from the 2D Fourier transform. 
Both approaches are performed in the present analysis to check for consistency. 

B. High-Speed Langmuir Probe Analysis 

High-speed Langmuir probe analysis is described in much greater detail in Lobbia’s work 25 and will only be 
briefly summarized here. The high-speed Langmuir probe data is first corrected for capacitive effects associated 
with high-speed sweeping. This is done by subtracting the null probe signal from the active probe signal. The data is 
then chopped into individual sweeps and analyzed using standard collisionless non-magnetized thin-sheath 
Langmuir probe theory. 26 Fourier transforms are performed on the resulting number density as a function of time to 
obtain power spectra. 


V. Results 


A. High-Speed Camera Results 

The discharge of the HiVHAc EDU2 was found to be dominated by breathing mode oscillations. For the 500 V 
throttle points at all test background pressures, there were no discernable signs of spokes mode oscillations. Spokes 
were found at other throttle points but were found to exist only within limited sectors around the thruster. 

Figures 6-9 show the theta-t diagrams (also called spokes surfaces) for the 500-3. 9-Pnx operating conditions. 
The theta-t diagram is a surface contour plot where the x-axis is the azimuthal axis and the y-axis is the time axis. 
Red indicates higher-than-average light intensity and blue indicates lower-than-average light intensity. Theta, the 
azimuthal coordinate, is defined as 0 degrees at the top of the annulus. Increasing theta is equivalent to clockwise 
movement on the image. Note the video is filmed through a mirror so the filmed image is a horizontally-mirrored 
image of the real thruster face and the reported theta coordinate is not corrected for this reflection. 



Theta, deg. 


Figure 6. Theta-t diagram for 500-3.9-Plx. 



Theta, deg. 

Figure 7. Theta-t diagram for 500-3.9-P3x. 
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Figure 8. Theta-t diagram for 500-3.9-P5x. 



Ql,,,- — . . 7 " ■ * 

0 45 90 135 180 225 270 315 360 

Theta, deg. 

Figure 9. Theta-t diagram for 500-3.9-P10x. 


From Figs. 6-9, one can see a large number of horizontal bands of color. These bands of color indicate a fairly 
uniform light intensity across all azimuthal locations that varies in time, which is the signature of breathing mode 
oscillations. In contrast, a spokes oscillation shows up as bands of diagonal lines because the light packet is traveling 
across both angle and time. Examination of 500-2.0 and 500-3.0 throttle point data as functions of background 
pressure indicates the same trend; there is a strong breathing mode of varying frequency present but no spokes mode 
oscillations. 

Figures 10-13 show the high-speed camera power spectra for the 500-3. 9-Pnx operating conditions with mode 
number 0 to 5 being shown. Note m = 0 corresponds to breathing mode while m = 1+ correspond to spokes mode 
with different modal patterns. The y-axis of the figures is spectral power on a logarithmic scale and x-axis is 
frequency in kilohertz. Since the videos were taken at 80,000 frames per second, the maximum resolved wave 
frequency is 40 kHz. 
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Frequency, kHz 

Figure 10. High-speed camera 
power spectra for 500-3.9-Plx. 



Frequency, kHz 

Figure 11. High-speed camera 
power spectra for 500-3.9-P3x. 



Frequency, kHz 

Figure 12. High-speed camera 
power spectra for 500-3.9-P5x. 
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Figure 13. High-speed camera 
power spectra for 500-3.9-P10x. 


While some peak structures can be seen for m = 1 and 2 the peak magnitude of these structures are two orders of 
magnitude lower than the m = 0 peak structures. The breathing mode oscillation is clearly dominant over spokes 
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mode oscillation at these operating conditions. Analysis of the 500-2.0 and 500-3.0 operating conditions show the 
same trend. 

Although high-speed videos for the background pressure test at the 300 V and 400 V throttle points are not 
available, there are videos at the lowest achievable pressure that were obtained during the performance testing phase 
of the component compatibility test. A search through these videos was performed to help determine why spokes 
mode oscillation was absent in the videos of the 500 V operation. Examination of thruster operation at 400-3. 2-Plx 
revealed a peculiar phenomenon where spokes mode is observed but only over a small section of the thruster 
annulus. Figures 14 and 15 show the theta-t diagrams and the high-speed camera power spectra, respectively, for the 
400-3. 2-P lx operating condition. 
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Figure 14. Theta-t diagram 
for 400-3.2-Plx. 


Figure 15. High-speed camera power spectra 
for 400-3.2-Plx. 


Figure 14 indicates that a single spoke repeatedly appears and disappears in the space between theta equal to 90 
and 135 degrees. The camera power spectra appear to indicate that the size of the spoke is a good match for the m = 
4 and 5 modal patterns. The frequency of the spokes is distinctly different from that of the breathing mode 
oscillation. Direct observation of the high-speed video confirms the presence of this peculiar spokes mode structure. 
Spokes mode structure does appear to form within the HiVHAc EDU2 but some damping mechanism damps out 
these spokes shortly after formation. Examination of the corresponding sector of the thruster post-test did not reveal 
any obvious differences compared to the rest of the thruster. While it is still unclear why spokes mode oscillations 
could not be sustained in the HiVHAc EDU2 operating under the tested conditions, it is clear that the analysis is 
capable of capturing spokes, and that spokes oscillations can be formed. The spokes oscillations are just not being 
sustained. 

Table 1 summarizes the breathing mode frequency and magnitude as functions of facility background pressure 
calculated from the high-speed videos. Oscillation magnitude is reported in the form of relative deviation in light 
intensity with a 95%-confidence-interval (95%-CI), much like how relative uncertainty is typically reported. The 
relative deviation is reported as a percentage of the time -averaged light intensity, where 95% of the time the light 
intensity falls within the reported boundary. For example, if 50% relative deviation is reported, that means 95% of 
the time, the light intensity fell within ±50% of the average light intensity. This table shows that the breathing mode 
frequency rises while the relative deviation of light intensity falls with rising background pressure. Recall from Eq. 
(1) the breathing mode frequency is theorized to scale with the square -root of the ion and neutral velocity and the 
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inverse of the characteristic length. At a given throttle point, retarding potential analyzer data previously reported 
shows that the ion velocity remains roughly constant as background pressure varies. 1 The neutral velocity scales 
with the square root of the anode temperature and is expected to remain roughly constant as background pressure 
varies. The implication is that the breathing mode frequency must be rising with background pressure because the 
characteristic length is decreasing. 

B. High-Speed Langmuir Probe Results 

Analysis of the high-speed Langmuir probe data focused on the breathing mode oscillations. Probe 
measurements taken at different axial locations are generally consistent in terms of measured oscillation frequency. 
Electron density and temperature are observed to decrease with axial distance while plasma potential asymptotes to 
low tens of volts. Only plasma data from the position 10 cm axially downstream of the exit plane will be reported 
since the whole data set is too large to show here. The important trends with background pressure can be observed 
with only the 10-cm-axial -position data. Data from both of the probe pairs at this position gave the same answer to 
within measurement uncertainty and the averaged values are presented. 

Table 2 summarizes the time-averaged electron density, electron temperature, and plasma potential at different 
operating conditions. The data are color-coded by pressure level for convenience. Generally speaking, the electron 
density increases with rising background pressure, which is expected if ingestion is occurring. Electron temperature 
and plasma potential both decrease with rising background pressure, suggesting a recession of the ion acceleration 
zone into the channel of the thruster. 

Table 3 summarizes the frequency of the plasma oscillation and the magnitudes of oscillations in electron 
density and electron temperature. Oscillation magnitude is reported in the form of relative deviation with a 95 %- 
confidence-interval. The oscillation frequencies of the electron density and electron temperature were found to be 
the same so only one frequency is reported. Some relative deviations of the electron temperature are greater than 
100%, which usually happens when there are periodic spikes in the electron temperature. This may have been an 
artifact of the data analysis routine or noise in the data. Despite this issue, the overall qualitative trends in the 
relative deviation of the electron temperature appear to be very consistent across throttle points. 

The oscillation frequency of the local plasma around the probe was found to be in fairly good agreement with 
that of the global plasma as measured by the high-speed video. This agreement gives us some certainty that the two 
diagnostics are both interrogating the same breathing mode phenomenon. The frequency calculated from the high- 


Table 2. Summary of time-averaged plasma properties measured by the high-speed Langmuir probe at 10 cm 

axially downstream of the exit plane along the thruster centerline. 

Disch. Disch. Electron density, m" 3 Electron temperature, Plasma potential, V 

voltage, power, eV 


V 

kW 

lx 

3x 

5x 

lOx 

lx 

3x 

5x 

lOx 

lx 

3x 

5x 

lOx 

300 

1.5 

6.4el6 

6.8el6 

- 

7.1el6 

10.1 

7.8 

- 

4.9 

53 

39 

- 

24 

300 

3.0 

1.3el7 

1.4el7 

- 

1.8el7 

5.5 

5.0 

- 

4.8 

23 

20 

- 

18 

400 

2.0 

5.3el6 

6.6e 16 

- 

7.1el6 

10.6 

9.5 

- 

6.0 

54 

43 

- 

27 

400 

3.2 

8.9el6 

9.4el6 

- 

1 . 1 el 7 

7.7 

6.4 

- 

4.5 

40 

32 

- 

21 

500 

2.0 

- 

3.4el6 

3.9el6 

- 

- 

7.4 

7.0 

- 

- 

36 

32 

- 

500 

3.0 

3.0el6 

5.4el6 

6.1el6 

7.2el6 

7.7 

6.3 

5.1 

3.7 

40 

33 

27 

20 

500 

3.9 

5.4el6 

7.2el6 

8.0el6 

8.9el6 

7.5 

5.6 

4.6 

3.5 

38 

30 

25 

18 


Table 3. Summary of time-resolved plasma properties measured by the high-speed Langmuir probe at 10 cm 
axially downstream of the exit plane along the thruster centerline. 


Disch. Disch. Relative deviation in Relative deviation in electron Oscillation frequency, kHz 

voltage, power, electron density, m 3 temperature, eV 


V 

kW 

lx 

3x 

5x 

lOx 

lx 

3x 

5x 

lOx 

lx 

3x 

5x 

lOx 

300 

1.5 

35% 

50% 

- 

61% 

26% 

58% 

- 

127% 

10.1 

10.7 

- 

11.2 

300 

3.0 

54% 

60% 

- 

81% 

183% 

227% 

- 

303% 

29 

37 

- 

40 

400 

2.0 

37% 

44% 

- 

63% 

31% 

44% 

- 

122% 

- 

33 

- 

38 

400 

3.2 

51% 

59% 

- 

64% 

56% 

95% 

- 

212% 

30 

- 

- 

39 

500 

2.0 

- 

32% 

35% 

- 

- 

29% 

37% 

- 

- 

25 

31 

- 

500 

3.0 

59% 

34% 

37% 

41% 

23% 

36% 

47% 

138% 

22 

28 

30 

31 

500 

3.9 

38% 

37% 

43% 

53% 

27% 

40% 

72% 

220% 

30 

29 

30 

36 
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speed video is of higher statistical significance because it was calculated from data around the entire annulus as 
opposed to a single location. Interestingly, the relative deviation in electron density and electron temperature 
generally increases with rising background pressure. This is the opposite trend as the relative deviation in light 
intensity from the high-speed video. 

One possible explanation is that the plasma downstream of the acceleration zone (where the Langmuir probes 
are) really does behave differently with changing background pressure than the plasma inside of the acceleration 
zone. Note that the high-speed camera signal is dominated by light emissions from regions of the plasma with high 
electron temperature, which is inside the ionization and acceleration zones of the thruster where electron 
temperature reaches tens of eV. Assuming the behavior of the plasma inside the channel differs from that of the 
plasma outside the channel, the combined results of the high-speed camera and Langmuir probes suggests a very 
interesting scenario. As the background pressure rises, the downstream boundary of the acceleration zone of the 
thruster is being pushed upstream, which is correlated with an increase in damping of the breathing mode magnitude 
inside the channel. At the same time, breathing mode magnitude outside the channel becomes less damped. If we 
further assume that the breathing mode oscillation enhances electron transport, the present data can be interpreted to 
mean that the shortening of the ionization and acceleration zones at high pressure reduces the amount of transport 
enhancement needed inside the channel (fewer magnetic field lines to transport across). At the same time, the 
increased distance between cathode and the downstream boundary of the acceleration zone at high pressure 
increases the transport enhancement needed to maintain the same flow of electrons to the channel. While this theory 
is consistent with the data, further experiments are needed to prove or disprove this theory. 

VI. Discussions 

A. Correlating Breathing Mode Frequency to Wear Mechanics and Wear Test Recommendations 

Conclusions about the change in the energy distribution of the ions bombarding the wall of the discharge channel 
as a function of background pressure cannot be made from the present study. However, one can correlate the 
breathing mode data to the length of the ionization/acceleration zone of the thruster and generate recommendations 
for testing from the results. 

As previously stated, assuming the scaling in Eq. (1) is correct, the breathing mode frequency only scales with 
ion velocity, neutral velocity, and characteristic length. For a given throttle point the ion and neutral velocities are 
fixed so the measured breathing mode frequency is an indirect measure of the characteristic length. Since the 
characteristic length is roughly the length of the ionization/acceleration zone, one can indirectly measure the relative 
change in the length of the ionization/acceleration zone as a function of background pressure. For a thruster like the 
HiVHAc EDU2, the part of the ionization/acceleration zone that is inside the discharge channel is roughly the same 
as the erosion zone. With a reasonable assumption about how the upstream edge of the ionization/acceleration zone 
is anchored, the breathing mode frequency can be used to determine how the upstream edge of the erosion band is 
shifting (the downstream edge is at the exit plane of the thruster). This information can be used to establish an upper 
limit on acceptable background pressure for which the change in erosion band length is considered acceptable. The 
same information can also be used to quantify the change in erosion band length. 

Figure 16 shows a plot of the normalized length scale as a function of normalized background pressure generated 
from high-speed camera data. Normalized length scale is the ratio of calculated length scale to the predicted length 
scale at zero background pressure. For example, a 
normalized length scale of 0.8 means that the calculated 
length scale is 80% as long as the length scale at zero 
background pressure. This figure shows a very consistent 
decrease in length scale as the background pressure is 
raised. 

To establish a background pressure limit for future test, 
a few assumptions will be made. For a typical Hall thruster, 
the acceleration/ionization zone length is a few tens of 
millimeters. ’ ’ Assuming a few millimeters change in 

the acceleration/ionization zone length correspond to a few 
millimeters change in the erosion band length and assuming 
the design margin on the erosion band length was chosen to 
be a few millimeters, one can establish 0.9 as the limit for 
how much change in normalized length scale is considered 
acceptable. Using Fig. 16, the normalized pressure 
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Figure 16. Normalized length scale calculated 
from high-speed video versus normalized 
background pressure. 


corresponding to a normalized length scale of 0.9 is about 2-3. Normalized pressure of 2.5 translates to about 1.2e-6 
Torr/(mg/s). In other words, the erosion band length will have changed by at most a few millimeters if the operating 
pressure is less than 1.2e-6 Torr/(mg/s). Note that this criterion does not account for changes in plasma energy and 
overall erosion rate, only erosion band length. For reference, the criterion on background pressure based on a 
maximum of 2% increase in thrust due to ingestion is 5e-7 Torr^mg/s). 1 

B. Updated Ingestion Model 

In a previous work, a two-stream ingestion model was 
proposed to account for ingestion of background neutrals in 
a Hall thruster. 1 The part of the model that predicts the 
properties of the external stream turned out to yield non- 
physical results leading to the conclusion that linearly 
super-imposing the internal (from propellant manifold) and 
external (from facility) streams on top of each other does 
not accurately model the ingestion physics. 1 Nevertheless, 
calculating the ingestion rate using the change in discharge 
current should be a valid approach since every ingested 
neutral will produce an electron upon ionization. The ion 
will be accelerated through a part of the potential difference 
between the anode and the cathode, and the electron will 
travel the other part of the potential on its way to the anode. 

Figure 17 shows an illustration of the two -stream ingestion 
model. For the present study, the model is modified so that 
the ingestion rate can be calculated without making any assumptions about the ingestion area. In the previous study, 
the ingestion area was assumed to be the same for all background pressures at each throttle point. 1 

The updated two-stream ingestion model still assumes two streams of neutrals. The first stream comes from the 
propellant manifold while the second stream comes from the facility. The streams have the same species 
composition, propellant utilization efficiency, and current utilization efficiency. Under these assumptions, the ratio 
of discharge current to mass flow rate for the ingested stream can be assumed to be equal to that of the main 
propellant stream. This ratio can be calculated using performance data from the component compatibility test. The 
discharge current is then plotted against the background pressure and a linear fit is performed. Zero-background- 
pressure performance of the thruster and ingested current can then be calculated. Ingested mass flow is then 
calculated from the ingested current using the aforementioned discharge -current-to-mass-flow-rate ratio. Table 4 
shows the measured thrust, discharge current, and calculated ingestion rate. Note that one entry in the calculated 
ingestion is negative. This unphysical data point is likely the result of a poor curve-fit or measurement uncertainty. 
Recall the accuracy of the mass flow controller is ±1% of reading. 


Va ^j KH* 

Ionization 

i i i 

◄- — Electron path 

Channel j 

Ion path 

Wall 



Anode flow 


Channel 

Wall 
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Downstream 
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Figure 17. Illustration of the simplified two 
stream ingestion model. V c 

voltages, not distances. 


and V ing indicate 


Table 4. Summary of ingestion rate calculation. 


Disch. Disch. Thrust, mN Discharge current, A Ingestion rate in percentage 

voltage, power, of supplied mass flow rate 


V 

kW 

lx 

3x 

5x 

lOx 

lx 

3x 

5x 

lOx 

lx 

3x 

5x 

lOx 

300 

1.5 

85 

90 

- 

95 

4.95 

5.06 

- 

5.31 

0.5% 

2.4% 

- 

6.9% 

300 

3.0 

188 

191 

- 

196 

10.03 

10.06 

- 

10.32 

0.5% 

0.8% 

- 

3.3% 

400 

2.0 

105 

109 

- 

110 

5.04 

5.26 

- 

5.24 

-0.9% 

3.1% 

- 

2.8% 

400 

3.2 

172 

175 

- 

181 

8.02 

8.05 

- 

8.19 

0.3% 

0.6% 

- 

2.3% 

500 

2.0 

92 

95 

96 

- 

4.06 

4.18 

4.19 

- 

0.4% 

3.1% 

3.3% 

- 

500 

3.0 

139 

144 

146 

149 

5.99 

6.08 

6.13 

6.22 

0.0% 

1.5% 

2.3% 

3.7% 

500 

3.9 

182 

186 

188 

193 

7.74 

7.77 

7.90 

7.96 

0.3% 

0.7% 

2.3% 

3.1% 


C. Dimensional Analysis on Rising Pressure and Shortening Ionization/Acceleration Zone 

While data from both the high-speed camera and the high-speed Langmuir probe suggests the 
ionization/acceleration zone is shortening and the downstream boundary is receding into the channel with rising 
background pressure, the root cause is unclear. A dimensional analysis was performed in an attempt to find the root 
cause. 
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Equation (2) shows the relationship between the plasma potential and charge density in the form of Poisson’s 
equation. In this equation, ® is the plasma potential, p is the charge density, and s 0 is the permittivity of free space. 
The corresponding scaling equation is shown in Eq. (3). Note that the Laplacian of ® scales as V/L 2 where V is 
voltage and L is length scale. 


V 2 O = -p/8 0 

T.~-L~ 1 _ 

>/MR a + MR j 

In Eq. (3), Li is the characteristic length scale, n e is the 
electron density, MR a is the anode mass flow rate, and MRi 
is the ingested mass flow rate. Figure 18 shows a plot of the 
normalized length scale predicted by the above dimensional 
scaling against ingestion percentage for the data set in the 
present study. This figure shows that the results under- 
predict the change in length scale by about one order of 
magnitude. A more complex approach is needed to 
accurately predict the change in transport that occurs with 
rising background pressure. 


VII. Conclusions 

The present study combined high-speed camera and 
high-speed Langmuir probe data to present a novel picture 
of the effect of facility background pressure on the plasma 
oscillations in the HiVHAc EDU2 Hall thruster. In the process of studying these oscillations, it was discovered that 
the discharge of the HiVHAc EDU2 at the operated conditions is dominated by breathing mode oscillations. Spokes 
can occasionally form but are quickly extinguished. 

The high-speed camera data showed that the breathing mode frequency varies from 10 to 36 kHz and tends to 
rise with rising background pressure. At the same time the magnitude of the light intensity oscillation inside the 
discharge channel falls with rising background pressure. Data from the high-speed Langmuir probe showed that the 
fluctuation in plasma density and temperature rises with background pressure 10 cm downstream of the exit plane. 
Data from the two diagnostics combined suggested that the ionization/acceleration zone of the thruster was 
shortening and receding into the discharge channel with rising background pressure. If true, the decreased 
characteristic plasma length would correspond to a reduction in the amount of field lines that electrons need to be 
transported across inside the discharge channel and may lead to a corresponding decrease in transport being 
provided by the breathing mode. This phenomenon would explain the decrease in the magnitude of light intensity 
oscillations. At the same time, the increase in distance between the downstream boundary of the acceleration zone 
and the cathode created an increased reliance on the breathing mode to carry the electron current from the cathode to 
the channel, which gave rise to an increase in fluctuation magnitude in plasma properties as seen by the Langmuir 
probes. 

The scaling relationship between breathing mode frequency and characteristic length scale was then used to 
establish a maximum acceptable background pressure for which the erosion band length will have changed by at 
most a few millimeters. This maximum acceptable pressure was found to be 1.2e-6 Torr per mg/s of anode mass 
flow rate. 

An update to the simple two-stream ingestion model was made to decrease its dependency on unknown 
quantities. Although the model is still not capable of accurately predicting the properties of the external stream, it 
can be used to calculate the rough ingestion rate without constraining the effective ingestion area. 

A simple dimensional analysis was performed regarding how the ingestion of background neutrals can cause the 
ionization/acceleration zone to shorten and recede into the discharge channel. Although the resulting prediction is 
qualitatively in the right direction, the prediction is off by roughly one order of magnitude. Further work is needed to 
determine the underlying mechanism that links rising background pressure with shortening ionization/acceleration 
zone. 

A number of assumptions were made in the course of analyzing the data presented in this paper. The most 
crucial assumption was that the scaling equation for the breathing mode frequency is correct. An important next step 



Ingestion Percentage 


Figure 18. Predictions of normalized length scale 
as a function of ingestion percentage. 


( 2 ) 

( 3 ) 
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will be to determine the validity of this scaling relationship. To experimentally verify the breathing mode frequency 
scaling relationship, one can deploy three diagnostics simultaneously. The first diagnostic will be the high-speed 
camera, which will be used to confirm the presence of the breathing mode and measure the associated frequency. 
The second diagnostic is laser-induced fluorescence velocimetry, which will be used determine the length and 
location of the acceleration/ionization zone as well as the velocity of the ions. The third diagnostic is either a set of 
thermocouples or a calibrated thermal camera to extract the anode temperature, which will be used to calculate the 
neutral velocity. The background pressure will then be varied at different discharge voltages and currents to produce 
the data needed to study the breathing mode frequency scaling relationship. 
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